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Abstract. We consider a simple stochastic differential equation for modeling bubbles 
in social context. A prime example is bubbles in asset pricing, but similar mechanisms 
may control a range of social phenomena driven by psychological factors (for example, 
O I- popularity of rock groups, or a number of students pursuing a given major). Our goal is 

to study the simplest possible model in which every term has a clear meaning and which 
demonstrates several key behaviors. The main factors that enter are tendency of mean 
reversion to a stable value, speculative social response triggered by trend following and 
random fluctuations. The interplay of these three forces may lead to bubble formation and 
collapse. Numerical simulations show that the equation has distinct regimes depending on 
• the values of the parameters. We perform rigorous analysis of the weakly random regime, 

' and study the role of change in fundamentals in igniting the bubble. 
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1. Introduction 

The best known and well-studied examples of social bubbles are speculative bubbles in 
asset pricing. These have a long history that in some form can be traced back at least to 
ancient Rome [TO]. First well-documented examples of speculative bubbles are the tulip 
q ■ mania in Netherlands in 1637 (see e.g. [131 EQl Elj) and the South Sea Company bubble of 

ON 

Generally, economists use the term "bubble" to describe an asset price that has risen above 
the level justified by the economy fundamentals, as measured by the discounted stream of 
expected future cash flows that would accrue to the owner of the asset. In practice, of 
course, what makes it difficult to diagnose a bubble is the uncertainty present in any 
^ \ calculation of the future returns and of the appropriate discount rates. 

Major modern examples of bubbles include the Japanese asset price bubble of 1980's [S] 
which involved both real estate and equities, the Dot.com bubble [9] in the US information 
technology stocks, and the recently punctured real estate and wider credit bubble [IB], 
which was centered in the US and UK but was present and affected a much wider range 
of markets. There have been many more smaller and local bubbles, involving regional real 
estate, stocks in a certain sector, and individual stocks [HJJ [29] . 

A simple intuitive explanation of the mechanism behind bubble formation has been sug- 
gested by Shiller [13]: "If asset prices start to rise strongly, the success of some investors 
attracts public attention that fuels the spread of the enthusiasm for the market. New 
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(often, less sophisticated) investors enter the market and bid up prices. This "irrational 
exuberance" heightens expectations of further price increases, as investors extrapolate re- 
cent price action far into the future. The market's meteoric rise is typically justified in 
the popular culture by some superficially plausible "new era" theory that validates the 
abandonment of traditional valuation metrics. But the bubble carries the seeds of its own 
destruction; if prices begin to sag, pessimism can take hold, causing some investors to exit 
the market. Downward price motion begets expectations of further downward motion, and 
so on, until the bottom is eventually reached" . 

It is reasonable to assume that bubble formation in asset prices has its root in some 
basic aspects of the human social psychology, which may manifest itself as soon as some 
necessary conditions (such as existence of basic liquid markets) are satisfied. Moreover, 
it is likely that bubble formation and collapse phenomena appear in areas of social dy- 
namics beyond the asset price formation. One example may be attraction of an excessive 
number of people to some profession, after there has been a significant breakthrough. The 
mechanism of the bubble in this case can be modeled very similarly to the price formation: 
explosive growth attracts newcomers who project growth into the future. Once the readily 
available new applications have been worked through, however, the profession may suffer 
from overcapacity of labor. Other relevant example can be the number of people following 
certain fashion or popularity of a rock group. 

The existing research literature on asset price bubble formation is too enormous to at- 
tempt a complete overview in this article. There are many different approaches and direc- 
tions. We only indicate a few branches. In finance literature, much work has been devoted 
to finite difference price evolution models. The first model of so-called rational bubbles 
based on finite differences was proposed by Blanchard and Watson [5]. The asset price sat- 
isfies a finite difference equation which basically expresses the no-arbitrage condition. The 
bubble component may be present due to the non-uniqueness of the solution. For discussion 
and more refined versions see also e.g. [TJ |2j [7J [191 EQl EH]. Continuous time models based 
on strict local martingale approach have been considered in [251 [2B1 EH E5] , where fur- 
ther references can be found. Another approach uses finitely additive measures (charges) 
to model the pricing bubbles (see e.g. [221 1221 I2Z])- The charges and local martingales 
approaches were recently shown to be equivalent in [28J. All the models discussed so far 
have no-arbitrage conditions explicitly built into them. 

Closest to our interest here is the direction of agent-based models of price formation. The 
primary motivation of this class of models is more realistic replication of the statistical prop- 
erties observable in the time series of stock prices returns. The agent-based models usually 
incorporate a number of agents following different strategies with effective price equation 
driven by the balance of supply and demand. The no-arbitrage conditions are typically 
absent, even though for sophisticated models constructing consistently winning strategies 
may be challenging. Two primary classes of strategies that received most attention are 
trend following ("chartist") and fundamentalist. The dynamics generated by interaction of 
trend following ("chartist") trading strategies and fundamentals-driven investors has been 
the focus of many works (see e.g. [U [HI EE CEO (2H1 [36]). Other authors pursued more 
sophisticated models based on heterogeneous and adaptive beliefs where strategies may 
vary in time (see e.g. [21 El E?l EHl EU H21 IZE]). Yet another direction of research is based 
on parallels with statistical physics and trading network models. An elegant theoretical 
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physics-inspired model leading to superexponential growth in prices while speculative heat 
lasts was proposed in [T5]. We refer to [HI US] for recent reviews of different directions and 
contributions, where many more references can be found. 

Our goal in this paper is to introduce and investigate a simple stochastic differential 
equation which models creation of bubbles in asset pricing. The equation has three key 
terms. The first is the mean reversion term, driving price back to the fundamental value. It 
models the contribution of the fundamentalist trading strategy, and is similar to what has 
been considered before. The second is random term, which models exogenous factors. The 
third term is the speculative or social response term, which models psychology of trend 
followers. Though the latter term is also similar in spirit to what has been considered 
before (see e.g. [Tl| [T7J EE]), the exact form is different in that our term is nonlinear, 
which is essential for certain dynamical features. 

The direction of our work is different from earlier literature in several respects. First, we 
intentionally keep the equation conceptually as simple as possible - though also sufficiently 
complex to produce rich set of phenomena. Thus we do not aim at this point to produce a 
model that can closely explain all of the observable statistical features of complex modern 
markets, but rather look for the simplest signature model of bubble and collapse, perhaps 
the next order of approximation to the reality after the random walk. One motivation is that 
even if not exhaustive, a simple model has a better chance of being capable of calibration 
to the empirical data without overfitting. There are only three independent parameters 
in the model, and we investigate the behavior of the model across the possible values of 
the parameters. In one of the regimes - small randomness - the model can be analyzed 
rigorously, and this provides valuable insight into the possible behavior in other regimes. 
Secondly, randomness plays more important role in our model than usual. In particular, in 
the absence of the change in fundamentals, randomness is entirely responsible for igniting 
the bubble and causing the bubble to collapse. This means that the deterministic part of 
our dynamics does not suggest any typical time scales for these processes, making them 
essentially random, and similar to Poisson process. Indeed, the bubble collapse (or ignition) 
are notoriously hard to predict. Third, we focus on the influence and role of the change in 
fundamentals in igniting the bubble, which is arguably a key reason behind the initiation 
of many bubbles. In our model's framework, a change in fundamentals may significantly 
increase the likelihood of bubble ignition. 

Here is the equation that we are going to study, written for the logarithm of the asset 
price: 

(1) dP(t) = -n(l - e p °- p V)dt + adB t + vS(P(t) - P(t - T))dt. 

The parameters /x, a and v regulate the strength of the mean reversion, random and 
speculative terms respectively, Pq is the fundamental price of the asset (that is assumed 
constant for the moment but will be allowed to vary later), B t is the Brownian motion 
and S(x) is the social response or speculative function. It will be assumed to have the 
following properties: S(x) is odd, monotone increasing and S(x) — > 1 as x — > oo. The 
function S has the natural structure reinforcing the existing trend in price dynamics. It is 
important, however, that the reinforcement strength depends on the rate of past returns 
in a certain way. We will assume that v ^> /i, so that for large values of the argument, 
the speculative function can dominate mean reversion. On the other hand, we will assume 
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that for small values of the argument, the speculative fever is negligible and is dominated 
by mean reversion; in particular, uS'(0) <C /i. 

The social response term of this structure looks reasonable for many problems where 
human psychology is concerned. In our opinion, it is reasonable that the size of the social 
response term varies nonlinearly depending on the strength of the trend. Indeed, weaker 
trends in pricing are not as eye catching and generate significantly less attention, news 
and press coverage. The exact shape of the function S, as far as it satisfied the properties 
outlined above, did not affect much the qualitative properties of the time series in our 
numerical experiments. One could argue, however, in favor of a more subtle dependence of 
the social response on the past trend. This would be equivalent to using a more complex 
memory integral operator as the argument for the social response function, rather than just 
P(t) — P(t — T). This more general case is interesting, and may be crucial for matching 
realistic price dynamics properties of liquid financial instruments (see the discussion section 
at the end of the paper for more details). However, in the current paper we restrict ourselves 
to a simple time delay for which equation ([!]) already exhibits rich behavior. 

Note that the parameters fx and v have dimension 1/t, and a has dimension l/yi. 
Without loss of generality, we can set the delay time T in (JTj) equal to one. Indeed, any 
time delay T can be reduced to this case by rescaling time and other coefficients. For the 
rest of this paper, we will fix T = 1. 

The model ([I]) has several essential regimes, in particular the stable mean reversion dy- 
namics, bubble, and collapse. These regimes are especially clear cut when the randomness 
is small, but become less evident when the randomness increases. The transition probabil- 
ities between the regimes depend on the parameters of the model. To develop intuition, in 
Section [2] we look at the possible regimes of the deterministic model with a = 0. The deter- 
ministic equation possesses the same key regimes as the random model, but no switching 
between the regimes is possible without randomness. In Section [3] we provide an heuristic 
simplified picture for the bubble equation behavior, which is useful to keep in mind when 
dealing with the general case. In Section HI we look at the random case. Here we show 
that, provided the randomness is small, our basic regimes remain stable with high proba- 
bility. We also derive estimates on the probability of switching between different regimes, 
which becomes possible with randomness. In Section we show that the probability of 
bubble creation can be enhanced greatly by manipulation of the stable value Pq. This may 
correspond to a strong earnings report exceeding expectations in the case of a stock, or 
to a stimulative interest rate policy in the case of the real estate market or commodities. 
In Section [BJ we present some basic numerical simulations, illustrating some of the results 
we prove and testing more general parameter values. A more extended set of simulations, 
in particular studying the statistical properties of the model (presence of fat tails for the 
distribution of returns? clustering of volatility? correlation functions?) are postponed 
to a later publication. In Section [7J we discuss various extensions and generalizations of 
our model addressing more realistic trend-following speculation term as well as modeling 
bubbles in spatially distributed systems. 
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2. The Deterministic Case and Regimes 

As a first step, we look at various regimes in the deterministic bubble equation dynamics, 
that is, we set a = 0. Let us introduce the shortcut notation 

/(P,P ) = -Ml-e P °- p(t) ). 

The properties of the function S will be very important, so let us recall once again that 
the function S is odd, increasing, and S(x) — > 1 as x — > oo. In addition to this, we will 
assume that there is no oscillation in S: more precisely, the second derivative of S(x) is 
continuous, positive on [0,6) for some b > and negative for x G (b, oo). The symmetry 
(oddness) assumption is not necessary and is made simply for convenience. The assumption 
of the saturation (existence of a finite limit of S(x) as x — > oo, rather than a continued 
growth for large values of x) may be debatable. Some authors suggest that often one can 
discern super-exponential growth approaching the height of a bubble [TTfl I2"5t 

Throughout the paper, we will make certain assumptions on the parameters and func- 
tions appearing in bubble equation ([1]). Our first assumption ensures this equation pos- 
sesses the mean reversion regime, where psychology has negligible effect on price dynamics. 

Assumption I. uS'(0) <C /i. 

We will need the following more technical version of this assumption to facilitate the 
proofs. 

Assumption P. There exists 5 m > such that we have 

f(P + x, P Q ) + vS(x + S m ) < -c m n5 m , for 5 m /4 < x < 5 m , 

and 

f(P + x, P ) + vS{x- 5 m ) > c m fx5 m , for -8 m /A >x> -5 m , 
where c m is a fixed positive constant. 

As we mentioned already in the introduction, Assumption I is quite natural. After all, if 
the trend is small, there is not much social excitement about it. It is reasonable to assume 
that the speculative fever starts only when the trend is significant. 

We now define the mean reversion regime and prove rigorously its stability in the absence 
of randomness. 

Mean reversion regime. Assume that \P{t) — Pq\ < 5 m for t G [to — lj^o]- Then the same 
holds true for any t > t . 

Proof. Assume that P(t) violates the corridor \P(t) — Pq\ < 8 m at some point and let ti 
be the minimal time greater than to when P(ti) — P$ = 5 m (the case P{ti) — Pq = —S m is 
similar). Then 

(2) P'(t) = f(P + 5 m , P ) + uS(P(t) - P(t - 1)) < f(P + 5 m , P Q ) + vS(28 m ) < 
by definition of t\ and Assumption I'. This is a contradiction. □ 

Figure [T] illustrates the mean reversion regime (the value of Pq is set to be zero here). 
The region between points £2 and £3 on this figure is the range of values of P where the 



(i 



ALEXANDER KISELEV AND LENYA RYZHIK 



Mean Reversion Regime 
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Figure 1. The graphs of f(P, 0) and —vS(2P). The region between X2 and 
£3 corresponds to the mean reversion regime. 

mean reversion term dominates the social response function, as in (j2J). The parameter 
5 m is the characteristic price scale of the mean reversion regime. Existence of such scale 
satisfying Assumption I' follows from Assumption I since df(Po,Po)/dP = — //. Below we 
will sometimes consider Pn depending on time. In that case, Assumption F holds at each 
t with the corresponding Po(t). 

Our second assumption ensures existence of the bubble and collapse regimes. 

Assumption II. Equation fi + x = uS(x) has exactly two positive roots, x 5 and x 6 , and 
exactly one negative root, x\. 

It follows that equation x = uS(x) has exactly two positive roots and exactly two negative 
roots, xi and £3 (due to symmetry). The smaller of these positive roots will be denoted 
x 4 . 

Let us explain why Assumption II is relevant for the presence of the bubble regime. 
Intuitively, a sustained bubble is driven by the balance between social response function, 
growth in P and the mean reversion term. For large values of P, the mean reversion term 
is basically equal to — //. Setting P' = x, and approximating P(t) — P{t — 1) m P'(t), we 
arrive at the balance 

(3) x + fi = vS(x) 

for the stable bubble regime. Given our assumptions on the structure of the response 
function S, equation (jSJ) can have two, one (this is a degenerate case) or zero positive 
roots. In the last two cases, the nonlinear response term is just too weak to sustain a 
bubble. Therefore, Assumption II is essentially just about the sufficient strength of the 
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Bubble and Collapse regimes: Case A 

50 r 




Figure 2. 

social response term. Observe that the smaller positive root, x$, is unstable, while x§ is 
stable (see Figure [2]). The root xq then determines the asymptotic rate of growth in a 
bubble. The root X5, on the other hand, determines the range of stability of the developed 
bubble regime (should a random fluctuation reduce the rate of growth below x$, the bubble 
will puncture and collapse will likely begin). Let us set 5 6 = (x$ — x 5 )/2 - this parameter 
can be regarded as the typical stability scale of the bubble regime (it has dimension of 
price/time). Note that it is clear from the definition that 25 > S m . 

The positive root £4 is relevant for values of P not too far away from Pq. When the mean 
reversion is small, X4 is the threshold rate of growth for the start of the bubble. However, if 
the rate of growth will not accelerate quickly to beyond x$, the beginning bubble may slow 
down and turn with the growth of P - due to the mean reversion term getting stronger 
and approaching the value —fi. Similarly to Xq, the stable negative root X\ determines the 
rate of collapse for large values of P. The collapse regime is extremely stable for large P. 
The system tends to go into a collapse with the rate of decay aproaching x\ if only the rate 
of growth falls below x 5 . 

Observe that without Assumption II, under our structure assumptions on the function S 
there is a possibility that (j3J) has three negative roots (or, in a degenerate case, two negative 
roots): see Figure [3] for illustration. In this case, the smaller (in absolute value) stable root 
Z2 determines the likely rate of decay right after the puncture of a bubble. We call such 
scenario "a non-panic collapse", since in many cases the primary driver of such decay is 
the mean reversion force that is close to the value — /x, and the social response function 
S does not play a big role. If \i is small, this regime may look like simply a pause in the 
bubble growth. A switch either into the bubble or into steeper collapse is possible if random 
fluctuations are present. Alternatively, the "non-panic collapse" can simply bring the price 
all the way down to the equilibrium value Pq, and switch to the mean reversion regime. 
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Bubble and Collapse regimes: Case B 




Figure 3. 

While this configuration is interesting and adds an extra phenomenon which deserves to 
be studied, in this paper we restrict ourselves to the simpler case covered by Assumption 
II: most collapses in practice tend to involve panic. The case with three negative roots will 
be studied elsewhere. 

For values of P closer to equilibrium, where f(P, Pq) cannot be regarded as just — /x, the 
stable negative root x 2 of % — vS(x) approximates the rate of collapse decay. The unstable 
root X3 provides a bound for the range of stability of the collapse (the collapse is definitely 
stable up to decay rates of 23; it really is more stable since the mean reversion works in 
favor of collapse for values P(t) > P ). We set 5 C = (x 3 —x 2 )/2, and regard 5 C as the typical 
stability scale of the collapse regime. Let a b = x 5 + 5b/ 2, and a c = x 3 — 5 c /2. Observe that 

(4) vS{a b ) - /i - a b = c b 5 b > 0, 

(5) uS(a c ) - a c = -c c 5 c < 0, 

where c b , c c > are some fixed constants depending on S. These inequalities, which follow 
from Assumption II, will be useful for us in the random case, where we will need some 
cushion to ensure the stability of the bubble and collapse regimes. 

We now describe the bubble and collapse regimes and their stability in the deterministic 
case. 

The bubble regime. Assume that P(t) — P(t — 1) > x§ for all t £ [to — l,^o]- Then we 
also have P(t) — P(t — 1) > 25 for all t > t . Moreover, P'{t) — > xq as t — > 00. 

Proof. Assume, on the contrary, that there exists a time t± such that P(t{) — P(t\ — 1) = X5, 
and t\ is the smallest such time greater than i - However in that case, for every t G [ti — 1, £1], 
we have 

P\t) > -fi + uS(P(t) - Pit - 1)) >-fi + uS(x 5 ) = x 5 



BUBBLE FORMATION AND COLLAPSE 9 

due to Assumption II. It follows that P{t\) — P(t 1 — 1) > x 5 , which is a contradiction. 
Hence, P'(t) > x 5 and 

(6) P(t) - P(t - 1) > x 5 for all t > t . 
Observe that, as a consequence, P(t) — > oo as t — )■ oo. Let us denote 

Pn = m i n t£[t +n-l,t +n]P{'t) > 

then the sequence P n is increasing by (J6]), and, moreover, P n ^ Co + 71x5, with some 
constant Co that depends only on the data on the time interval [to — l,to]- We also define 

p+ = max te[to+n _i ito+n ]P'(i), 

and 

p~ = mm te[t0+n _ lit0+n] P'(t). 

We claim that 

(7) lim p~ = lim p\ = x 6 . 

n— >oo n— >oo 

To this end we make the following observations. First, note that p~ > x$ and p^ > x§ for 
all n. We also claim that 

(8) if p~ < x 6 then p~ +1 ^ p~ . 

Indeed, assume that p~ ^ x 6 . For any t E [t + n, t + n + 1] we have 

(9) P'(t) = -/i + M e p °- p + uS{P{t) - P(t - 1)) > -/i + ^(min(p;,p; +1 )), 
thus 

(10) p- +1 > -n + isS(mm(p-,p- +1 )). 
Now, if (ED fails then implies that 

(11) p~ +1 > -fi + uS{p- +1 ). 

As p n+ \ > x 5 it follows from (TlTl) that > x 6 , hence we have the inequalities 

which is a contradiction. Thus, ([8]) holds. 
Next, we show that 

(12) if p^ ^ Xq then p~ ^ x 6 for all n ^ k. 

It suffices to show that (Tl2l) holds for n = k + 1. However, it follows from ( flOl) that either 

(13) > -/i + ^(p* ) ^ + uS(x 6 ) = x 6 , 
or 

( 14 ) Pk+i > -V + v S (Pk+i)> 
which also implies that p^ +1 > x 6 . Therefore, f|T2l holds. 

The next step is to show that there exists a constant cq so that 

(15) if p+ ^ x 6 + c e" nx ' 5 then p+ +1 ^ p+. 
Note that 

pt +1 < -/i + jxe flj - p " + i/S(max(p+ «C -/i + /ie^^e"™ 5 + */£(max(p+ p++i)), 
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hence 

(16) p+ +1 <-fi + d e- nX5 + isS{max{p+,p+ +1 )), 

with do = /xexp(P — Co)- Now, if p+ < Pn+i then 

p+ +1 <-^ + d e- nX5 + vS(p+ +1 ), 

which, in turn, implies that 

vt < Pn+l < ^6 + C e- nX5 

with some constant Cq that depends only on the behavior of the function S(x) near x = Xq, 



and ( 1151) follows. 

Now, let us prove (GO). Let us set 

L + = lim sup 

Passing to the limit in (I16p gives 

L + ^ -n+ vS{L + ), 
hence L + ^ xq. Then ([8]) and (Tl2|) imply that 

L~ = lim p~ 

n— >oo 

exists and £5 < L~ ^ x 6 . Moreover, (flOj) implies that 

L" ^ -fi + u3(L-), 

whence L~ = x 6 . As L + $C x 6 , and L + ^ it follows that L + = x 6 as well, and, 
moreover, 

lim inf p+ is lim inf p~ = Xq , 

n—^oo n— >oo 

thus © holds. □ 

While the mean reversion and bubble regimes in the deterministic case continue indef- 
initely, the collapse regime ends once we reach the equilibrium value (in fact, as will be 
clear from the argument below, we should go well below the equilibrium value to break the 
decay) . 

The collapse regime. Assume that both P(t) — P(t — 1) < x 3 and P(t) > P for all 
t G [to — 1, to]- Then for any t > to such that P(s) remains larger than Pq for all to < s < t, 
we have P(t) - P(t - 1) < x 3 . Moreover, if P(t) > P for t e [to - 1, t + 1], and 

x 3 > maxt 6[t0 _i it0 ](P(t) - P(t - 1)) = d > x 2 , 

then the decay accelerates in the next time interval: 

maxt e[ i oA5+ i](P(t) - P(t - 1)) = di< d . 



Remark. We lack a statement on the asymptotic rate of collapse since any collapse ends 
once Po is reached. But if one assumes that P is very large in the beginning, one can prove 
that the rate of decay approaches X\ while P remains very large. The argument is similar 
to that given above in the bubble case. 
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Proof. Suppose there exists ti > t such that for t G [to,ti), we have P(t) — P(t — 1) < x 3 
and P(t) > P , but P(t x ) - P(ti - 1) = x 3 . Then for any t G [h - l,ti], we have 

P'(t) < uS(P(t) - P(t - 1)) < uS(x 3 ) = x 3 , 

giving us a contradiction. 

Now let max tg [ to _i jto ](P(i) — P(t — 1)) = do, with X2 < d < x 3 . Then for every t G 
[to — 1, to], we have P'(t) < vS(do) < d by the definition of X2 and x 3 . This implies, by a 
familiar argument, that P(t) — P(t — 1) < d for every t G [to, to + 1]- D 

It is clear that P should fall sufficiently below P to balance the speculative function. 
In particular, if the fall has been from large values of P, we need at least /(P, Po) > x 2 to 
break the fall. 

We should stress that it is clear from our conditions that mean reversion, bubble and 
collapse do not exhaust the whole range of possible behaviors in our model. For most 
parameter ranges where the system does exhibit interesting behavior, there is a significant 
gap between mean reversion and bubble. We call it a transitory regime. In the random 
case, one situation where it appears is when a fluctuation kicks the price away from P 
but not strongly enough for the bubble to form (for example, the rate of growth does 
not exceed x^ > 5 m or is just barely above it). In this regime, there is no dominating 
term. The effects of both mean reversion and social response may be significant. In the 
deterministic setting, what happens depends strongly both on the rate of growth and the 
value of P: for example, if P > Po and max t6 [ to _ lito ](P(t) — P(t — 1)) < x±, the growth 
will decelerate. In the random case, the evolution in the transitory regime becomes much 
more sensitive to randomness. The transitory regime is thus most difficult for rigorous 
analysis. In this paper, we will prove rigorous bounds on probabilities reflecting stability 
or transitions between mean reversion, bubble and collapse. The bounds we prove may be 
not quite sharp due to the transitory effects we just described. 

Observe that in the deterministic case, the three basic regimes are not connected: one 
cannot transition between them, with the exception of the collapse regime which ends once 
the price has fallen enough. Other transitions will be made possible by randomness. There 
is a way, however, to ignite bubbles in the deterministic case as well - by adjusting the level 
of the equilibrium value, P . This can be interpreted as a change in the fundamental data: 
a new positive earnings report in the case of a stock or a stimulative interest rate policy for 
real estate or fixed income investments. Since in practice, such change in P often happens 
in a jump, that is what we will consider here. It is not difficult to adjust the arguments 
below to the case of a continuous dependence Po(t). 

The first observation is that if the jump in P is small, then the bubble does not ignite. 

Driving deterministic bubbles: a small jump. Assume that |P(t) — Po(t)| < S m /2 for 
t G [to — l,to], Po{t) — P-- Assume also that Po jumps up at to to a value P + and the 
size of the jump satisfies < P + — P_ < S m /2. Then we have P(t) — Po(t) < S m /2 and 
P(t) — Po(t) > — 5 m for t G [to, t + 1], Po(t) = P+. Thus we continue in the mean reversion 
regime for all future times. 

Remark. We need a more relaxed lower bound above simply because the original bound 
may clearly fail right after the jump. 
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Proof. Assume on the contrary that P(ti) = P + + 5 m /2 at some t 1 for the first time. Then 
we have 

P'ih) = f(P + + SJ2, P + ) + vSiPfa) - P(h - 1)), 

and, moreover, 

P(f - P(h - 1) < P^) -p + + p + -p_ + p_- P( tl - l) < 5 m /2 + 5 m /2 + 5 m /2 = 3<W2. 

Using Assumption F with x = S m /2, gives 

/(P + + 8J2, P + ) + uS(35 m /2) < 0, 

and so P'(ti) < which is a contradiction. Similarly, assume that Pfa) = P+ — S m for 
some t2 > to for the first time. But then 

P'{h) = f{P + - 5 m , P + ) + uS(P + -S m - P(t 2 - 1)) 

> f(P+ - Sm, P+) + vS(P+ -5 m -P_- SJ2) > 

by Assumption F. □ 

On the other hand, a strong driving will necessarily lead to bubble ignition. 

Driving deterministic bubbles: a large jump. Assume that \P(t) — Po(t)\ < S m /2 for 
t E [t — l,t ], Po(t) = P— Assume that P jumps up at t to a value P + and the size of 
the jump satisfies 

P+ - P- > x 5 + — + 5 m . 

H 

Then for all t E [t + 1, t + 2] we have P(t) - P(t - 1) > x 5 . 

Proof. First, observe that for all t G [t , t + 1] we have P(t) — P(t — 1) > — 5 m . This is true 
since P(t — 1) < P_ + S m /2 by assumption, and P(t) > P — S m /2 by the same argument 
as in the proof of stability of mean reversion regime (using P + > P_). 

Now we claim that P'(t) > x 5 , for every t G [to, io + !]■ Indeed, if P(t) < P_ + x$ + 5 m /2, 
then 

P'{t) > f(P(t), P + ) + uS(-6 m ) = Me p +- p W - 1) - Ke Sm/2 - 1). 
We used Assumption F in the last step. Then 

P\t) > MP+ - P(t) - 5 m /2) > x 5 . 

On the other hand, if P(t) > P_ + x 5 + 5 rn /2, then 

P'(t) >-fi + uS(P(t) - P_ - 5 m /2) >-fjL + vS(x 5 ) > x 5 . 

Now we can prove that P(t) — P{t — 1) > x 5 for all t G [to + 1, *o + 2] by a familiar 
argument. Indeed, this is true for t = t Q . Assuming existence of some minimal t\ where we 
have P{ti) — P(ti — 1) = x 5 , we quickly obtain a contradiction. □ 
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3. The random case: heuristics 

In this section, we will provide an intuitive picture of what to expect from our model (0Q) 
when the random term is included, but is small. This regime may be more relevant to 
modeling illiquid investments like real estate rather than liquid ones. Consider first the 
mean reversion regime. In this case we can regard the influence of the social response term 
as minor. It gets activated only once randomness, by accident, moves the price away from 
the mean reversion zone. Note that if we gradually drift away from Pq - which is difficult 
due to mean reversion - we still do not activate the bubble regime as the social response 
term is small, but move into the transitory regime. Rather, we should have a jump of size 
~ a b to enter the stable bubble regime. Let us set P n = P(t + n), AP n = P n — P„_i, 
AB n = B n — B n _i. Then heuristically, the evolution of P n in mean reversion regime can 
be approximated by 

AP n = aAB n -(t{P n -P ). 
Intuitively, we pass into the stable bubble regime on the nth step if aAB n > (1 + fi)at,. 
Thus the probability Ai of bubble ignition at each step is roug hl y 1 - $(H±hK) ) w here $ 
is the cumulative function of the normal distribution, 

$(ar) = -)= f e- y2/2 dy. 

V27T J-oo 

At this heuristic level, the dynamics in the mean reversion regime can be thought of as a 
Poisson-type process where on each time step we pass to the bubble regime with probability 
Ai, and stay in the mean reversion regime with probability 1 — Ai- Of course, in reality 
the probability of bubble ignition depends on the past dynamics and on what exactly 
happens within the time interval. The social response term also increasingly plays a role 
as the system moves into the transitory regime which is more prone to creating bubbles 
than mean reversion. However, the above picture captures the most qualitative aspects of 
dynamics. 

Once we are in the bubble regime, our rate of growth approaches xq, and the bubble 
stability in terms of variations in the rate of growth is measured by 5^. The dynamics is 
qualitatively described by 

(17) AP n = -n + a AB n + uS{AP n ^). 

Now, in order for the bubble to burst and for dynamics to pass to the collapse regime, 
we need aAB n > 5b- Thus, the dynamics in the bubble regime can be modeled by a 
Poisson-type process where on every time step we switch from bubble to collapse with 
the probability A 2 ~ 1 — Q(5b/a) and we stay in the bubble with the probability 1 — A 2 . 
Observe the similarity between this heuristic picture and the well known simple discrete 
model proposed and investigated by Blanchard and Watson [5] (minus the no-arbitrage 
condition, which is not automatic in our model). 

Finally, once we have entered the collapse regime, the dynamics in the leading order 
is described again by (1171) (for larger values of P, with the term (— u) on the right side 
gradually disappearing as P decreases). For large P, the only way to break the collapse 
is to go back into the bubble regime, and for that one needs uAB n > a& — a c . Observe 
that the probability of this is small compared to Ai and A 2 , due to a& > 5f,, a c < —5 C . 
As we approach the equilibrium value Pq, the condition for breaking the collapse gradually 
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approaches aAB n > S c , and the passage to the transitory or mean reversion regime becomes 
possible. Once we have broken through P , the push back to the mean generated by 
the mean reversion term increases rapidly We pass from collapse to the recovery when 
P -P~log(l + |). 

We emphasize that the heuristic picture of this section cannot be expected to provide 
a close approximation of the real evolution in all possible relevant ranges of parameter 
values. For example, when the random term is sufficiently strong, transitory effects play 
an important role, and dependence of the evolution on details of the past dynamics increases 
notably. This regime is perhaps most interesting from the practical point of view. Still, in 
Section E] devoted to numerical simulations, we will see that our heuristic picture appears 
to describe the evolution pretty well for a significant range of parameter values. 

4. The Random Case: Mean Reversion, Bubble and Collapse 

Now, we consider the full equation ([I]): 

dP = f(P)dt + adB t + uS(P(t) - P(t - l))dt. 

Naturally, depending on the relative strength of parameters, this equation can exhibit very 
different behavior. For example, if both v and a are small, the mean reversion term will 
always dominate and we will have a process very close to Ornstein-Uhlenbeck process. Tak- 
ing u large emphasizes randomness, and then the nonlinear social response effects may be 
difficult to discern. In this section we provide rigorous analysis which is most relevant in 
the case of small to moderate random forcing. First, we describe the probabilistic stability 
of our basic regimes. For now, the equilibrium value Pq is fixed. 

Stability of the mean reversion regime. Assume that for all t G [to ~ l,to] we have 
\P(t) - P | < $m, and \P(t ) - P | < <W2. Then with probability at least 2$(C5 m /a) - 1, 
we have that the same is true for t G [to,t + 1] : \P(t) — Po| < S m for t G [to,t + 1] and 
\P(to + 1) — Po| < S m /2. Here C > is a universal constant. 

Remark. A stronger condition at the end of the interval is necessary: if we had |P(to) — Pol — 
5 m , then we would have exited the interval [Po — 5 m , Pq + 5 m ] at later times with probability 
one. 

Proof. Let us assume that for any t\ < t2, ti, £2 G [to? to + 1] we have 

(18) a\B t . z - B tl \ < min(c m /i, l/4)5 m , 

where c m is the same as in Assumption I'. This is true in particular if 

crmax to < t < to+1 |P t - B to \ < ^min(/ic m , 1/A)6 m . 

The probability of this event is equal to 2$(^ 2i ) — 1. 

We claim that in this case, |P(t) — Po| < 5 m for all t G [t , t + 1] and \P(to + 1) — Po| < 
5 m /2. Let us verify that P(t) < P + 5 m for all t G [t , t + 1] and P(t + 1) < P + 5 m /2. 
The other part of the condition can be verified similarly. Indeed, assume on the contrary 
that there exists t 3 with P(t 3 ) = P + 5 m . Let t 4 be the largest time in [to,ti] such that 
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P(t 4 ) = P + S m /2. Then we have 

ta 

P(t 3 ) - P{U) = J (f(P(s), P ) + uS(P(s) - P(s - 1))) ds + a(B t3 - B u ). 

t4 

Note that due to P(s — 1) > Po — 5 m and Assumption F, we have 

f(P(s), P ) + uS(P(s) - P(s - 1)) < -c mf i5 m 

for any s in the interval of integration. Therefore we get 

P(t 3 ) - P(t 4 ) < -c m /i<5 m (t 3 - U) + 5 m /A < 5 m /4, 

a contradiction to our choice of £3, £4. 

Similarly, under condition ( ITS]) we have P(t + 1) < P + S m /2. Indeed, assume that 
P(t + 1) > Po + <W 2 - If we ever have P(t) < Po + for t e [t , t + 1], let t 3 be the 
largest time less than t + 1 when P^) = Po + <5 m /4. Then 

P(t + 1) - P(t 3 ) = J (f(P(s), P ) + ^(P(s) - P(s - 1))) rfs + «r(B t0+1 - S t3 ) 



< -c m /i5 m (t + 1 - h) + -r < ^ 

which is a contradiction. We used Assumption I' in estimating the integral. On the other 
hand, if P(t) > P + 5 m /4 for all t G [to 5*0 + 1] (but, as we saw above, P{t) < P + S m ), 
then 

t +l 

P(t + 1) - P(t ) = f (f(P(s), Po) + uS(P(s) - P(s - 1))) ds + (7(P t0+1 - B t0 ) 



to 

*^ c m /i5 m -|- c m <5 m /i 0, 

once again a contradiction. 

Thus, under our condition on the maximum of Brownian motion, the mean reversion 
regime is preserved. □ 

Stability of the bubble regime. Assume that for t G [to — 1, to], we have P(t) — P(t—1) > a&, 
and (7maxt 6 [i _i >to ]|.B i — -B to _i| < Q,<5(,/4. Then with probability at least 2$(C<5f>/cr) — 1, the 
same is true for t G [to,t + 1] (but to — 1 is replaced by t ). 

Proof. Suppose that for any t\ < t 2 , ti,t 2 G [to, to + 1] we have cr(B t , 2 — B tl ) > —CbSb/2. 
This is true with probability ta least 2$(C<5b/cr) — 1. Let us show that in this case, 

Pit) - P{t - 1) > a b for all t G [t ,t + 1]. 

This holds for t by our assumption. Let ti be the first time when P(ti) — P(ti — 1) = 
Then 

Pit,) - P(t x - 1) > -pi + vSia b ) + a(B tl - B tl _ x ) > a b 
under the assumption on Brownian motion and (j4]). This is a contradiction. □ 
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Remark. We need an extra assumption on the behavior of Brownian motion during 
the past interval [t — 1, to] since there are some unlikely scenarios where, due to Brownian 
motion, we essentially set up the exit from the bubble regime during [to — 1, to] even though 
P(t) — P(t — 1) stays large for this time interval. 

Stability of the Collapse regime. Assume that P{t) — P(t — 1) < a c and P(t) > Po for 
t G [t — l,to]. Assume also that amax t e[t -i,t ]\(Bt — Bt -i\ > c c 5 c /A. Then with probability 
at least 2$(C5 c /er) — 1 either there exists a time t G [to, to + 1] so that P(t) < Po ,or 
P(t) - P(t - 1) < a c for all t G [t , t + l}. 

Proof. Let us assume that for any t\ < t 2 , ti,t 2 G [t , t + 1] we have o~(B t2 — B tl ) < c c 5 c /2. 
This is true with probability exceeding 2Q(C5 c /a) — 1. Assume, on the contrary, that 
P(t) > P for all t G [t , t + 1], but there exists ti such that P(t x ) - P{ti - 1) = a c . But 
then 

P(h) - P(h - 1) < vS{a c ) + a{B tl - B tl ^) <a c - c c 5 c + c c 5 c /2 < a c , 
which is a contradiction. □ 

As in the previous case, it is easy to construct a scenario showing that we need to as- 
sume something about the Brownian motion on the interval [to — l,to] i n order to obtain 
a reasonable bound on the probability of the continuation of the collapse regime. 

Next, we look at a new phenomenon which appears in our model due to randomness. 
The bubble ignition probability. Assume that for t G [to — 1, to], we have \P(t) — Po\ < 8 m , 
and |-P(to) — -Pol < Then with a positive probability exceeding 

(19) po = e -C^ + * m+c ^)>' ^ ^ min(cUc^) X _ ^ ^ ^ _ ^ 
we have P(t) - P(t - 1) > a b for all t G [t + 1, t + 2]. Here c b is as in (jU). 



Proof. We will identify just one scenario leading to bubble ignition, and will estimate its 
probability. Assume that 

(20) \a(B t - B t0 ) - (a b + 2^ + 5 m + c b 5 b )(t - t )| < ^< 6 ^c b 5 b ) ^ 

for all t G [to, t + 1] (one sample path satisfying this bound appears on Figure H]). 

We will estimate the probability of ( 120]) in Lemma [1] below. First, we claim that under 
assumption ( 120]) . we have P(t) — P(t — 1) > — 25 m for all t G [to, to + 1]. This follows 
from P(t — 1) < P + 5 m due to our assumptions, and from P(t) > P — 5 m . The latter 
inequality can be established, given (120]) . similarly to the proof of the stability of random 
mean reversion regime. 

Secondly, we claim that P(t + 1) — -P(to) > a b- Indeed, 

P(t + 1) - P(t ) > -/x + a{B t0+1 - B t0 ) + vS{-25 m ). 

Observe that due to Assumption I', uS(— 25 m ) = —vS{25 m ) > —ji. Then using (|2TJ|) . we 
obtain 

P(t + 1) - P(t ) > -2/i + (a b + 2fi + 5 m + c b 6 b ) - min (^» c ^) > a& . 
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Figure 4. One possible behavior of the Brownian motion leading to the 
bubble ignition 

Finally, assume that -P(ti) — P(t\ — l) = a b for the first time at some t\ G [t + 1, to + 2]. 
Let us split 

P(h) - P(h - 1) = P{t x ) - P(to + 1) + P(to + 1) - P(h - 1). 

Then 

P(h) - P(t + 1) > -//fa - t - 1) + - B t0+l ) + ^(a 6 )(t! - t - 1) 

> (a b + c b 5 b )(ti - t - 1) + cr(S tl - B to+1 ). 

We used (TjJ in the second step. Assume that for any t 2 < t 3 , t 2 , t 3 G [t + 1, t + 2] we have 
a(B t3 — B t2 ) > —c b 5 b /4. This happens with probability exceeding 2$ ( £ ^ L ) — 1. Then 

(21) P(h) - P(t + 1) > (o 6 + c b 5 b ){h - t - 1) - oA/4. 

Next, 

to+l 

P(t + l)-P(t 1 -l)>- l i(t + 2-t 1 ) + a(B to+1 -B tl ^) + u J S{P{t)-P{t-l))dt 

ti-i 

> -2n(t + 2 - tx) + (2/i + a 6 + 5 m + qA)(t + 2 - t x ) - mm ^' > 

(22) > (a b + 5 m + c b 5 b )(t + 2 - t x ) - min ^' Cfe ^ . 

Adding up ([21) and ([22), we obtain 

P(tl) - P(ti - 1) > Oft + oA/4 > a 6 , 
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and this is a contradiction. 

It remains to estimate the probability of (|2"U1) . 

Lemma 1. Let X t = B t — ct. Then 

(23) P(max < s < t |X,| < b) > e- cb ^{2^(b/Vi) - 1). 

Proof. The proof is a simple application of the Girsanov theorem. Let A be the set of 

c 2 t 

paths where max <s<t|^ s | < b. Let Z t = e cBt ~ , and define a new measure dQ t = Z t dP t . 
Then by Girsanov's theorem, X t is a Brownian motion with respect to dQ t , hence Q(A) > 
2$(6/y / t) — 1. Denote by \a the characteristic function of this set, and by Eq the expec- 
tation taken with respect to the measure dQ. Then 

P{A) = E Q [e- cB ^ XA ] > E Q [e-< ct+b ^ X A} = e~ cb -^Q(A). 
We used that on A, B t < ct + b. Thus, 

P(A) > e- cb ~^r (2$(b/y/i) - 1 



□ 



It follows from Lemma [T] that the probability of (1201) satisfies 
P I max < s <i 



a b + 2/i + 5 m + c b 5 b 
B.« 1 



min(^ m , c b 5 b ) > 



e -(a b +2[i+6 m +c b S b ) 2 /<T 



2 ^ ^ mm{5 m ,c b 5 b ) ^ _ ^ 



(observe that the cb term in the exponent in f )23|) is dominated by c 2 t/2 in our case). 
Multiplying it by the probability of an independent Brownian motion behavior we required 
on [t + 1, t + 2] leads to (USD- □ 

Next, we look at the probability of the bubble regime switching to a collapse. 
The bubble collapse probability. Assume that we are in a Bubble regime: for t G [to - 1, ^o], 
P(t) — P(t — 1) > a b and P(t) > Po- Then with probability exceeding 

(24) e-C-- +^)> 2 (2* (|^) - l) (2$ ( min( ^ C ^^ - 

we have either P(t) < P for some t e [t ,t + 2], or P(t) - P(t — 1) < a c for all t e 
[t + l,to + 2]. 

Proof. The argument is similar to the bubble ignition case, and we leave it to the interested 
reader to complete it. □ 

One of the debatable features of the bubble equation is that the bubble to collapse 
transition probability does not appear to increase with growth of P. In fact, there is a 
slight increase due to decay in f(P,P ), but it is likely that the effect is insignificant 
when v ^> \i. In the discussion section we propose a slightly more complex model that 
is much more sensitive to the level of P. The drawback is the presence of an additional 
parameter controlling the negative feedback response to growth in price. Given how hard 
this parameter is to estimate in practice, we believe there is some merit in keeping things 
simple and elegant as in the bubble equation, where the timing of collapse is purely random. 
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5. The Random Case: Forcing a Bubble 
In this section, we consider the random case with a changing Pq. 

Driving random bubbles: a small jump. Assume that for all t G [to — 1 , to] ? we have 
\P(t)-P (t)\ < S m , \P(t )-P (t)\ < 5 m /2 where P (t) = P_. Assume that at t the stable 
value P undergoes a jump up to P + , with P + — P < <5 m /4, and Po(t) stays equal to P + 
after that. Then with probability at least 2$(C5 m /o") — 1, we have \P(t) — P+\ < 5 m for 
t G [t , t + 1], and |P(t + 1) - P+\ < S m /2. 

Proof. Let us estimate the probability of P(t) > P + — 5 m for t G [t , t + 1] and P(t + 1) > 
P + — 5 m /2. The other half of the statement is similar and is in fact a bit easier since there 
is more danger in violating the lower boundary due to the increase in Pq. Similarly to 
the stability argument without change in P , assume that o-\B t . A — B t2 \ < min(c m /i, l/4)S m 
for any t 3 > t 2 , t 2 ,t 3 G [to, to + 1]- The probability of this exceeds 2$ (|^) — 1. Assume 
that there exists t\ such that P(ti) — P + = —S m , and look at the minimal such t±. Since 
P(to) — P+ > — 3<5 m /4, find the largest t 4 G [to,ti] such that P(t 4 ) — P + = — 3£ m /4. Consider 

(25) P(tx) - P(t 4 ) = J /(P(t), P + ) alt + v J S{P{t) - P(t - 1)) alt + a(B tl - B u ). 
Observe that 

P(t) - P(t - 1) > P(t) - P + + P + - P_ - 5 m /2 > P(t) - P + - 35 m /4. 
Then by Assumption I', we have 

/(P(t), P + ) + z/S(P(t) - P(t - 1)) > c mf i5 m 
in ( 1251) . Given our Brownian motion assumption, we obtain 

P(ti) - P(t 4 ) > -5 m /4 + c my u5 m (ti - t 4 ) > -8 m /A, 

a contradiction. 

Next, assume that P(t + 1) < P + — 5 m /2. If we ever have P(t) > P + — <5 m /4 in [t , t + 1], 
let t 5 be the largest time where we have P(t 5 ) = P + — <5 m /4. Then 

/•to+l 

P(t + 1) - P(t 5 ) = / (f(P(s), P+) + ^(P(s) - P(s - 1))) rfs + 
Jt 5 

cr(B to+ i - B tb ) > c m fi5 m (t + 1 - t 5 ) - 8 m /A, 

giving a contradiction. We used Assumption I' in the second step. In the case if P(t) < 
P + — <5 m /4 (but we saw P(t) > P + — 5 m ) for all t G [to, t + 1], we have 

P(t + 1) - P(t ) > c m /i<5 m - c m /jS m > 0, 

contradiction. □ 

Now we consider the probability of bubble ignition given a strong jump in P . 
Driving random bubbles: a large jump. Assume that for all times t G [t — 1, to] we have 
|P(t) - P (t)| < 5 m , and |P(t ) - P (t )| < 6 m /2, with P (t) = P_. Assume that at t = t 
Po jumps to a value P + , and the size of jump satisfies 

(26) P + - P_ > 2Qfe + 3 ^ . 

min(/^, 1) 
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Then with probability exceeding 2$(Cmin(5 m , c b 5 b )/a) — 1, we have P(t) — P(t — 1) > a b 
for alltG [t + 1, to + 2]. 

Proof. The main step in the proof is the following claim. Suppose, in addition to the 
assumptions above, that for any t 2 ,t 3 E [to,t + 1], t 3 > t 2 , we have 

(27) a{B t3 - B t2 ) > -min(5 m , c b S b )/2. 
Then 

(28) P{h)-P{t 2 ) > (a b + mm(5 m ,c b 5 b ))(t 3 -t 2 )-min(5 m ,c b 5 b )/2. 

To prove (1281) . first observe that under our assumptions P(t) — P(t — 1) > — 25 m for all 
t E [to, to + 1]- This can be verified using (12 7p similarly to the proof of stability in random 
mean reversion regime. Now 

(29) P(t 3 ) - P(t 2 ) = f 3 (/(P(t), P + ) + vS{P{t) - Pit - 1))) dt + a(B t3 - B t2 ). 

If, on the other hand, P(t) — P_ > a b + S m , then P(t) — P(t — 1) > a b and therefore 
uS(P(t) - P(t - 1)) > a 6 + qA + /i. Since /(P(t), P+) > -/i for all values of P(t), we find 
that the expression under integral in (|29|) is greater then Ofc+Cb^. Now if P(t)— P_ < af>+<5 m , 
then by (|26|) we have 

P+ _ P(t ) > 

min(//, 1) 

In this case the expression under integral in f[2^|) exceeds 

^ e srR^T) _ i + vS(-26 m ) > fi e=wi - e a ™ > /i f 7 " > a b + 5 m . 
V / V / mm(/i, 1) 

Here in the first step we used uS(— 25 m ) > /i(l — e 5m ), which follows from Assumption I'. 
Combining these estimates and assumption fl27|) leads to f )28|) . 

Next, suppose that the Brownian motion on the interval [to + 1, to + 2] satisfies 

(30) a(B t3 - B t2 ) > -c b 6 b /2 

for any t 3 ,t 2 E [t + l,t + 2], t 3 > t 2 . Observe that implies that P(t + 1) - P(t ) > a b . 
Let us show that the inequality persists for all t E [t + l,t + 2] given (I3"U|) . Assume, on 
the contrary, that ti G [to + 1, to + 2] is the minimal time such that P(ti) — P(ti — 1) = a b . 
Split P(ti) - P(ti - 1) = P(tx) - P(t + 1) + P(t + 1) - P(ti - 1). Then 

P(t x )-P(t +1) > (vS(a b )-/j l )(t l -to-l)+o-(B tl -B to+1 ) > (a 6 +C6(y 6 )(ti-to-l)-C6V2. 

On the other hand, (128|) implies that 

P(t + 1) - P(ti - 1) > (a b + min(5 m , c b 5 b ))(t + 2 - t x ) - min(5 m , c b 5 b )/2. 

Combining these two estimates, we find that P(ti) — P(ti — 1) > a b - contradiction. 
It remains to observe that the probability of both f[2T|) and fl3"U|) happening can be 

estimated from below by 2$ ( min ^^ b) ) - 1- □ 
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6. Numerical Simulations 

We have tested the bubble equation dynamics on a wide range of parameter values. 
The Matlab random number generator was used for simulations. The model seems to be 
quite robust to the type of random forcing: mean zero uniformly distributed and normally 
distributed random terms with comparable variances lead to very similar results. This is 
again to be expected from the analytical point of view, since the properties of the Brownian 
motion used in the proofs are fairly generic and shared by many other random processes. 
The simulations were performed with the nonlinear response term S(x) = arctan(cir 2n+1 ), 
with n = 1, 2 or 3 and d a varying parameter. The roots Xi, i = 1, ...6 were computed for 
each parameter set tried. When clear bubble and collapse regimes were observable, the 
rates of growth and decay showed very good agreement with the ones predicted by values 
of the relevant roots. 

As one can expect, the behavior of the model depends strongly on the balance between 
different parameters. First of all, if we want to observe bubbles, we should set v > \i. We 
found three distinct regimes for the Bubble equation. 




Figure 5. Big rare bubbles phase: > aa > cr 2 

Regime 1: v > \i > a 2 , or big rare bubbles regime. The three regimes described in 
analytical part of the paper are very clear cut here. Transitions are relatively rare. The 
heuristic picture of Section [3] seems to apply very well. A typical simulation is shown on 
Figure |5j In this regime, we noted that the likelihood of serial bubbles occurence, passing 
from complete collapse right to the next bubble, may be significant. This happens if the 

price compression below Pq (which according to the heuristics is of the order ~ log 11 + —]) 
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Figure 6. Big rare bubbles phase: v > \i > a 2 : Serial bubbles. 

is comparable to x$. In that case the recovery to Pq level provides enough of a trend to 
help ignite the next bubble (see Figure [6]). The run on Figure |5] corresponds to v = 5, 
H = 4, a = 3, S cubic with d = 0.4. This corresponds to X\ ~ —12, x 2 ~ —7.5, x 3 ~ —0.5, 
£4 ~ 0.5, £5 ~ 2, x 6 ~ 3, and <5 m ~ 1. It is clear from these values that the stability of the 
mean reversion and bubble regimes is actually even stronger than one can expect from our 
analytical bounds. In fact, taking smaller a usually led to very long mean reversion runs, 
followed by eventual bubbles which grew huge and exceeded our computer capacity before 
collapsing. 

The fact that the dynamics in big rare bubbles regime is described quite well by heuristics 
of Section [3] makes it also close to the well known bubble model of Blanchard and Watson 
[5]. It would be of interest to show that under certain scaling assumptions, the dynamics 
of the bubble equation converges in the rigorous sense to a Poisson-type process with just 
three regimes of constant, exponential growth and exponential collapse. 
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Regime 2: a 2 > v > /i, or strong randomness regime. In this phase the effects of 
the social response term are weak compared to random fluctuations. The mean reversion, 
bubble and collapse regimes are difficult to isolate, and the dynamics is far from the heuristic 
model of Section |3j A sample simulation of this regime is shown on Figure [7J One graph is 
the price evolution given by the bubble equation, and the other, provided for comparison, 
is the dynamics corresponding to v — (and so it is just an Ornstein-Uhlenbeck process). 

The bubble equation price tends to overshoot Ornstein-Uhlenbeck, but fairly slightly, 
and there is no sustained bubble regime. The parameters are v = 0.6, fi = 0.2, o = 20, the 
function S is quintic with d = 90. 



14 




Figure 8. The balanced phase: a ~ v > fi 

Regime 3. The most interesting phase appears to be when a is slightly larger or compa- 
rable to u, and fi is also comparable or weaker. It appears that this is the parameter range 
that is most likely to be relevant in most applications with liquid instruments. Here one can 
still sometimes distinguish the different regimes (although mean reversion is very unstable 
if n is small). Yet transitory effects are stronger than in big rare bubbles phase. Figures |8] 
and M show some typical runs in this regime. As before, the Ornstein-Uhlenbeck process 
corresponding to the same random sample but with v = is graphed for comparison. 

Figure |8] corresponds to v = 0.6, fi = 0.23 and a = 2. The function 5* is quintic, with 
d = 90. The Figure M corresponds to v = 0.42, fi = 0.15, a = 2, S is cubic with d = 90. 

Figure M corresponds to relatively stronger randomness. While the nonlinear effects are 
still strong, fairly long transitory regimes become possible. Figure El shows an interesting 
phenomenon of a collapse regime switching back into bubble significantly above Pq- It is 
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Figure 9. The balanced phase: stronger randomness 

clear by looking at the corresponding Ornstein-Uhlenbeck process graph that randomness 
drives these switches. In both simulations, social response term increases the effective 
volatility of dynamics, stronger on the upside (creating bubbles) but also on the downside 
after collapse (this difference in the strength of the effect has clearly to do with very quick 
growth of the mean reversion function once the price is below Pq). 

Next we investigate the likelihood of a jump in fundamentals igniting a bubble. We 
find that the probability of the bubble ignition becomes significant if the jump becomes 
comparable in size to v. The Figure [TD] shows a simulation where the jump does not cause 
a bubble, while Figure HT1 shows a larger jump causing a bubble to ignite. On both figures, 
P(t), Po(t) and the Ornstein-Uhlenbeck process corresponding to v = case are graphed. 

The final graph shows a simulation where the stable value Po depends on time (Figure [12]) . 
The particular model we used is a random walk with a drift; the volatility of Po is much 
lower than that of P. All other parameters of the bubble equation are the same as in the 
run corresponding to Figure [9) The graphs on Figure [T2l are those of Po{t), Pit) and the 
corresponding Ornstein-Uhlenbeck process. 

7. Discussion and Further Directions 

The main goal of this paper was to present a simple differential equations model of the 
effects of human psychology in asset pricing. We tried to keep the model as simple as 
possible conditioned on the fact that it should exhibit sufficiently rich set of behaviors. A 
very interesting question that we did not address in this paper is the statistical properties 
of price returns generated by ([!]). It is well known that empirical daily price returns of 
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FIGURE 10. A small jump in fundamentals: no bubble ignition 

various financial instruments such as for example major indexes demonstrate deviations 
from normal distribution, in particular power-like tails in the distribution of the returns. 
Typically, lack of strong correlations in the ACF of daily returns is observed (even though 
for some indexes and for some historical periods significant correlation can be present 
(see e.g. |33j). On the other hand, the volatility of returns exhibits clustering with ACF 
decaying at a slow power rate. A number of models have been proposed that can fairly 
accurately account for these statistical properties (see e.g. [37J HSl HI] , and also JUJH^ for 
reviews). These models, however, tend to be more sophisticated than the one considered 
here. Can one obtain similar results for a simpler model ([I])? It is likely that some 
modifications will be needed, in particular to reduce correlations in the return series. We 
plan to discuss the statistical properties of the model and some of its modifications in a 
future publication. 

Another possible future direction of research involves bubble formation in spatially dis- 
tributed systems, such as the real estate markets. The price is now a function of both 
space and time, P(x,t). In the simplest case, the problem can be set on a graph, where 
each vertex Vj, j = 1,...,N corresponds to a town or a city. The interaction between 
prices pi and pj in cities Vi and Vj is defined by a diffusion coefficient Ky. This coefficient 
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Figure 11. A large jump in fundamentals causes bubble ignition 

models the degree of contact these cities have, geographical proximity and generally the 
degree to which the real estate prices in one city are likely to influence the prices in the 
other. To each Vi we also assign a number q^, which measures the size of v^s market. In 
addition, the stable value p° varies from city to city. Then one possible model for price 
dynamics is given by 

N 1 

dPi(t) = Ki i( pi ~ + f(Pi'P°i) + adB t + vSiPiit) - Pii 1 - !))■ 
3=1 aj 

The relative simplicity of the equation ([1]) describing price in every graph vertex (city) 
makes the model look approachable. Some of the interesting questions to study in this 
case are possibility and likelihood of bubble contagion (or bubble front propagation) doe 
to price diffusion, or possibility of existence of localized bubbles. 

Yet another interesting problem is suggested by numerical simulations in the big rare 
bubbles regime. Is it true that under certain scaling assumptions, the equation ([I]) can 
be rigorously shown to converge to a nontrivial limit, such as suggested in Section [HP 
Does one get simply three possible clear cut regimes with Poisson-like switching between 
them? A result like that would establish an analytical link between relatively sophisticated 
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FIGURE 12. The balanced phase: varying Pq 

bubble equation and simpler discrete bubble models reminiscent of the one appearing in 
Blanchard- Watson [5]. 
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